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ABSTRACT 

We proposed a new method - Pacman - to calculate Faraday rotation measure (RM) maps 
from multi-frequency polarisation angle data (Dolag et al.) in order to avoid the so-called nn- 
ambiguity. Here, we apply our Pacman algorithm to two polarisation data sets of extended 
radio sources in the Abell 2255 and the Hydra A cluster, and compare the RM maps obtained 
using Pacman to RM maps obtained employing already existing methods. Thereby, we pro- 
vide a new high quality RM map of the Hydra north lobe which is in a good agreement with 
the existing one but find significant differences in the case of the south lobe of Hydra A. We 
demonstrate the reliability and the robustness of Pacman. In order to study the influence of 
map making artefacts, which are imprinted by wrong solutions to the nTt-ambiguities, and 
of the error treatment of the data, we calculated and compared magnetic field power spectra 
fr om various RM maps. The power spectra were derived using the method recently proposed 
by EnBlin & Vogt ( 2003). We demonstrate the sensitivity of statistical analysis to artefacts and 
noise in the RM maps and thus, we demonstrate the importance of an unambiguous determina- 
tion of RM maps and an understanding of the nature of the noise in the data. We introduce and 
perform statistical tests to estimate the quality of the derived RM maps, which demonstrate 
the quality improvements due to Pacman . 

Key words: Intergalactic medium - Galaxies: cluster: general 



1 INTRODUCTION 

Magnetic fields are ubiquitous throughout the universe. One 
method to study them is the analysis of the Faraday rotation ef- 
fect. This effect describes that if polarised radiation traverses a 
magnetised plasma the plane of polarisation of the radiation is 
rotated. The change in polarisation angle is proportional to the 
squared wavelength A 2 . The proportionality constant is called rota- 
tion measure (RM). Various methods have been proposed in order 
to ca lculate RM maps from multi-frequency polarisation data sets 
(e.g .lVallee & Kronbergll975l : lHaveJl975l:lRuzmaikin & SokolofJ 
i 2001). The difficulty involved in this calcula- 



Il9y*isaraia& Jain 
tion arises from the fact that observations constrain the polarisation 
angle ip only up to additions of ±nn. This leads to the so-called 
nn- ambiguity. 

hereafter Paper I), we 



In a previous paper ( Dolag et al. 20 
proposed a new approach for the unambiguous determination of 
RM and of the intrinsic polarisation angle ip° at the observed, po- 
larised radio source from multi-frequency polarisation data sets. 
The proposed algorithm uses a global scheme in order to solve 
the nn- ambiguity. It assumes that if regions exhibit small polari- 
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sation angle gradients between neighbouring pixels in all observed 
frequencies simultaneously then these pixels can be considered as 
connected. Information about one pixel can be used for neighbour- 
ing ones and especially the solution to the ?m-ambiguity should be 
the same. 

Faraday rotation measure maps are analysed in order to get in- 
sight into the properties of the RM producing magnetic fields such 
as field strengths and correlation lengths. Artefacts in RM maps 
which result from nn- ambiguities can le ad to misinterpretation of 
the data. Recently, EnBl in & Vogtl fc003h proposed a method to 
calculate magnetic power spectra from RM maps and to estimate 
magnetic field prope rties. They suc cessfully applied their method 
to RM maps (Vogt & EnBlin 2003) and realised that map making 
artefacts and small scale pixel noise have a tremendous influence 
on the shape of the magnetic power spectra on large Fourier scales 
- small real space scales. Thus, for the application of these kind 
of statistical methods it is desirable to produce unambiguous RM 
maps with as little noise as possible. This and similar applications 
are our motivations for designing Pacman (Paper I) and to quantify 
its performance (this Paper). 

In order to d etect and to estimate the correlated noise level in 
RM and <p° maps, EnBlin et alJJ2003h recently proposed a gradient 
vector product statistic. It compares RM gradients and intrinsic po- 
larisation angle gradients and aims to detect correlated fluctuations 
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on small scales, since RM and ip° are both derived from the same 
set of polarisation angle maps. 

In Paper I, we describe in detail the idea and the implemen- 
tation of the Pacman 1 algorithm. Furthermore, we presented a test 
of the algorithm on artificially generated maps where we demon- 
strated that Pacman yields the right solution to the nir- ambiguity. 
Here, we want to apply this algorithm to polarisation observa- 
tion data sets o f two extended pola rised radio sources located in 
the Abell 2255 (iGovoni et all hoOD cluster and the Hydra cluster 
dTavlor & Perlevlll993h . This is described in Sect.0 Using these 
polarisation data, we demonstrate the stability of our Pacman al- 
gorithm and compare RM maps obtained from a standard fit al- 
gorithm. A "standard fit" algorithm is considered to be the algo- 
rithm which performs a pixel by pixel lea s t squa r es fit u s ing th e 
methods suggested by IVallee & Kronberd (Il975h : iHaveJ Jl975h : 
|Ruzmaikin&SokoMilil979|). 

We demonstrate the importance of the unambiguous determi- 
nation of RMs for a st atistical analysis by ap plying the statistical 
approach developed by Enfili n & VogJ fc003l) to the RM maps in 
order to derive the power spectra and strength of the magnetic fields 
in the intra-cluster medium. We also discuss the influence of error 
treatment in the analysis. The philosophy of this statistical analy- 
sis and the calculation of the power spectra is briefly outlined in 
Sect. l2.3l whereas the results of the application to the RM maps are 
presented in Sect. [3] 

In additio n, we apply the gra dient vector product statistic V 
as proposed by Enfil in et alJ fc003h in order to detect map making 
artefacts and correlated noise. The concept of this statistic is briefly 
explained in Sect. 12. H and in Sect. [3] it is applied to the data. After 
the discussion of our results, we give our conclusions and lessons 
learned from the data during the course of this work in Sect. [4] 

Throughout the rest of the paper, we assume a Hubble constant 
of H = 70 km s" 1 Mpc -1 , tt m = 0.3 and tt A = 0.7 in a flat 
universe. The notation used follows Paper I. 



2 STATISTICAL ANALYSIS OF RM MAPS 
2.1 Gradient Vector Product Statistic 



Enfilin et al. I fc003h introduced a gradient vector product statistic V 
to reveal correlated noise in the data. Observed RM and ip° maps 
will always have some correlated fluctuations on small scales, since 
they are both calculated from the same set of polarisation angle 
maps, leading to correlated fluctuation in both maps. The noise 
correlation between if and RM errors is an anti-correlation of lin- 
ear shape and can be detected by comparing the gradients of these 
quantities. A suitable quantity to detect correlated noise between 
if and RM is therefore 

/ d 2 xVRM(x) ■ Vif°(x) 



V 



J d 2 x\VRM(x)\\Vip°(x)\ 1 



(1) 



where VRM(x) and V(^°(cc) are the gradients of RM and ip°. A 
map pair, which was constructed from a set of independent random 
polarisation angle maps <^(/c), will give sl V > —1. A map pair 
without any correlated noise will give a V ~ 0. Hence, the statistic 
V is suitable to detect especially correlated small scale pixel noise. 

The denominator in Eq. (lj is the normalisation and enables 
the comparison of V for RM maps of different sources, since V is 

1 The computer code for Pacman will be publicly available at 

|http : //dipastro .pd. astro . it/ ~cosmo/Pacman| 



proportional to the fraction of gradients which are artefacts. Since, 
we especially want to compare the quality between maps of the 
same source calculated using the two different algorithms, it is use- 
ful to introduce the unnormalised quantity 



V = J d 2 xVRM(x) • V(p°(x) 



(2) 



The quantity V gives the absolute measure of correlation between 
the gradient alignments of RM and ip°. The smaller this value is the 
smaller is the total level of correlated noise in the maps. 



2.2 Error Under or Over Estimation 

One problem, we were faced with during the course of our statisti- 
cal analysis is the possibility that the measurement errors are under 
or overestimated. Such a hypothesis can be tested by performing a 
reduced xl test which is considered to be a measure for the good- 
ness of each least squares fit and calculates as follows 
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where v — f — n c is the number of degrees of freedom and / is the 
number of frequencies used and n c is the number of model param- 
eters in our case, n c = 2. The parameter s 2 is the variance of the 
fit and (crfc. ,)k is the weighted average of the individual variances: 



(crkij)k — 



(4) 



If one believes in the assumption of Gaussian noise, that the 
data are not corrupted and that the linear fit is an appropriate model 
then any statistical deviation from unity of the map average xl of 
this value indicates an under or over error estimation. For a xl > 1, 
the errors have been underestimated and for a xl < 1> they have 
been overestimated. 

Unfortunately this analysis does not reveal in its simple form 
at which frequency the errors might be under or overestimated. 
However, one can test for the influence of single frequencies by 
leaving out the appropriate frequency for the calculation of the RM 
map and comparing the resulting x^.. values with the original one. 

It is advisable to evaluate the x^.. -maps in order to locate re- 
gions of too high or too low values for xl by eye. 



2.3 Magnetic Power Spectra 

Our Pacman algorithm is especially useful for the determination of 
RM maps of extended radio sources. RM maps of extended radio 
sources are frequently analysed in terms of correlation length and 
the RM producing magnetic field strength. One analysis relying on 
statistical methods in order to derive magn e tic fie ld power spectra 
was recently developed bv lEnBlin & Vogtl y003). They success- 
fully applied their method to existing RM map s of radio sources lo - 
cated in Abell 2634, Abell 400 and Hydra A iVogt & EnfiliiJ2003h . 
In the course of their work, they realised that their method is sen- 
sitive to map making artefacts and pixel noise, which can lead to 
misinterpretation of the data at hand. Therefore this method is a 
good opportunity to study the influences of noise and artifacts on 
the resulting power spectra. 

The observational nature of the data is taken into account 
by introducing a window function, which can be interpreted as 
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the sampling volume. Especially, a noise reducing data weight- 
ing scheme can be introduced to account for observational noise. 
One reasonable choice of weighting is to introduce a factor 
a^ M /a^ M , which we call simple window weighting. Another pos- 
sible weighting scheme is a thresholding scheme described by 
1/(1 + a^ M /a^ M ), which means that the noise below a certain 
threshold a^ M is acceptable and areas of higher noise are down- 
weighted. 

Note, that magnetic power spectra, which are calculated using 
the approach explained above, are used as valuable estimator for 
magnetic field strengths and correlation lengths. However, since the 
spectra are shaped by and are very sensitive to small scale noise and 
map making artefacts, we use them here to quantify the influence 
of noise and artefacts rather than as an estimator for characteristic 
properties of the RM producing magnetic fields. 



3 APPLICATION TO DATA 

For the calculation of what we call further on the standard fi t RM 
maps, we used the y 2 m ethod suggested by IVallee & Kron berg 
Jl975h and lHaveJ Jl975h . Also for the individual RM fits, we 
adopted this method. We used as an upper limit for the RM values 
\BM\ij = RM max ± 3crg M , where crg M was calculated follow- 
ing Eq.(5) in Paper I. The least- squares fit for the individual data 
points were always performed as error weighted fits as described 
by Eq. (3) and (4) in Paper I, if not stated otherwise. For the Pac- 
man calculations, we used for the free parameters a = ft = 1 in 
Eq. (7) in Paper I. 

However, we a lso used the method suggested by 

iRuzmaikin & Sokolof i |l979) as standard fit algorithm. The 
results did not change substantially. 

3.1 Abell 2255E 

The Abell cluster 2 255, which has a redshift of 0.0806 
JStruble & Roodll999h . has been studied bv lBurns et aTNl995h and 
Fere tti et al.Nl997h . The polarised radio source B 1713+641, which 
we call hereafter A2255E, has a two sided radio lobe structure and 
is not directly located in the cluster centre. Polarisation observa- 
tions were performed using the Very Large Array (VLA) at 4535, 
4885, 8085, and 8465 MHz. The data reduction was done with stan- 
dard AIPS (Astronomical Imaging Processing Software) routines 
(Govoni et al., in prep.). The polarisation angle maps and their stan- 
dard error maps for the four frequencies were kindly provided to us 
by Federica Govoni. A preliminary analysis of an RM map of this 
source in order to determine the properties of the magnetic field 
in the intra-cluster gas in Abell 2255 is presented in lGovoni et alJ 
J2002h . 

Using the polarisation angle data, an RM map employing the 
standard fit algorithm was calculated, where the maximal allowed 
error in polarisation angle was chosen to be cr™ ax = 25°. The re- 
sulting RM map is shown in the upper left panel of Fig [Q The 
overlaid contour indicates the area which would be covered if po- 
larisation angle errors were limited by cr™ ax = 10°. For these cal- 
culations, we assumed RM max = 1 500 rad m~ 2 . 

As can be seen from this map, typical RM values range be- 
tween -100 rad m -2 and +210 rad m -2 . Sinc e this cluster is no t 
known to inhabit a cooling core in its centre iFeretti et alJll997h . 
these are expected values. However, the occurrence of RM values 
around 1000 rad m -2 , which can be seen as grey areas in the stan- 
dard fit map of A2255E, might indicate that for these areas the 



n7r-ambiguity was not properly solved. One reason for this sus- 
picion is found in the rapid change of RM values occurring be- 
tween 1 or 2 pixel from 100 to 1000 rad m -2 . All these jumps 
have a ARM of about 1000 rad m -2 indicating ?m-ambiguities 
between 4 GHz and 8 GHz which can be theoretically calculated 
from ARM = 7t/(A| — A 2 ). Another important point to note is 
that parts of the grey areas lie well within the cr™ ax = 10° con- 
tour, and hence, contribute to the results of any statistical analysis 
of such a parametrised RM map. 

Therefore, the polarisation data of this source provide a good 
possibility to demonstrate the robustness and the reliability of our 
proposed algorithm. The RM map calculated by Pacman is shown 
in the middle left panel. The same numbers for the corresponding 
parameters were used; a^ ax = 25°, RM max = 1 500 rad m~ 2 . 

An initial comparison by eye reveals that the method used in 
the standard fit procedure produces spurious RMs in noisy regions 
(manifested as grey areas in the upper left panel of Fig. [I) as pre- 
viously mentioned while the RM map calculated by Pacman shows 
no such grey areas. The apparent jumps of about 1000 rad m -2 
observed in the standard fit map have disappeared in the Pacman 
map. Note, that in principle by realising that these steps are due to 
n7r-ambiguities, one could re-run the standard fit algorithm with a 
lower RM m ax and this wrong solutions would disappear. However, 
by doing this a very strong bias is introduced since there might also 
be large RM values which are real. This bias can be easily relaxed 
by using Pacman. 

A pixel-by-pixel comparison is shown in Fig. In this fig- 
ure the RMstdfit values obtained for pixels using a standard fit are 
plotted on the y-axis against the RMp aC man values obtained at the 
same corresponding pixel locations using the Pacman algorithm on 
the x-axis. The black points represent the error weighted standard 
fit whereas the red points represent the non-error weighted standard 
fit. Note that Pacman applied only an error weighted least- squares 
fit to the data. From this scatter plot, one can clearly see that both 
methods yield same RM values for most of the points which is 
expected from the visual comparison of both maps. However, the 
points for RM s tdfit values at around ± 1000 rad m -2 are due to 
the wrong solution of the nn- ambiguity found by the standard fit 
algorithm. 

For the demonstration of the effect of nn- ambiguity arte- 
facts and pixel noise on statistical analysis, we determined the 
cluster magnetic field power spectra by employing the approach 
briefly described in Sect. 12.31 For detailed discussion and descrip- 
ti on of the application o f the method to data, we refer the reader 
to IVogt & EnBlinl fc003h . In the calculation, we assumed that the 
source plane is parallel to the observer plane. Please note, that the 
application of an RM data filter in order to remove bad data or data 
which might suffer from a wrong s olution to the nn- ambiguity as 
described bv IVogt & EnB lin (2003) is here not necessary and not 
desirable. This is on the one hand due to the use of an error weight- 
ing scheme in the window function to suppress bad data and on the 
other hand due to the wish to study the influence of noise and map 
making artefacts on the power spectra. 

The RM area used for the calculation of any power spectrum 
is the RM area which would be covered by the Pacman algorithm 
while using the same set of parameter cr™ ax , RM max . This en- 
sures that pixels on the border or noisy regions of the image not 
associated with the source - as seen in the standard fit RM map of 
A2255E in the upper left panel of Fig. [I]- are not considered for 
the analysis. (The same philosophy also applies for the discussion 
of the Hydra data in Sect. l3.2l and Sect. 13.31 . 

The power spectra calculated for various map making see- 
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Figure 1. Comparison of the standard fit RM maps (upper panels) and the Pacman RM maps (middle panels) for A2255E (left column), the north (middle 
column) and the south lobe (right column) of Hydra A. RM values not represented by the colour bar are coloured in grey. Panels at the bottom exhibit the 
respective cr RM -maps. For the standard fit RM map of A2255E a cr™ ax = 25° and k m i n = f = 4 was used, contours indicate the area covered by a standard 
fits RM map with cr™ ax = 10°. The parameters used for the Pacman RM map of A2255E are cr™ ax = 25°, cr^ ax = 25°, g = 1.2 and k m i n = f = 4; 
for the standard fit RM maps of Hydra A are cr™ ax = 25° and k min = f = 4 (at 8 GHz); for the Pacman RM map of Hydra (north) are cr™ ax = 30°, 
= 4, f = 5, a* = 25° and a = 1.5; for the Pacman RM map of Hydra (south) are cr™ ax = 35°, k m . im = f = 5, <r* = 35° and a = 2.0. 



Pacman II: Application 5 



A2255E 

2000 f 




-2000 I ] 

-1000 -500 500 1000 

^Pocmon 

Figure 2. A pixel-by-pixel scatter plot of RM stc jfit values obtained using 
the standard fit on the y-axis versus the RMp acman values calculated from 
the polarisation data of Abell 2255 by employing our new algorithm Pac- 
man. Black points represent results from an error weighted standard least- 
squares fit and red points from a non-error weighted standard fit. The points 
for RMgtdfit values of around ± 1000 radm -2 are artefacts of a wrong 
solution of the ?m-ambiguity. 

narios are shown in Fig. [5] The solid and dashed line represent 
power spectra calculated from standard fit RM maps. The solid 
line was calculated from an RM map which was obtained by using 
cr™ ax = 20° where no error weighting was applied to the standard 
fit. The possible window weighting as introduced in Sect. 12.31 was 
also not applied. This scenario is therefore considered as a worst 
case. The dashed line was calculated from an RM map employing 
an error weighted standard fit while only allowing errors in the po- 
larisation angle of cr™ ax = 10°. Again no window weighting was 
applied. From the comparison of these two power spectra alone, 
one can see that the large k- scales - and thus small real space scales 
- are sensitive to pixel-noise. Therefore, the application of an error 
weighted least square fit in the RM calculation leads to a reduction 
of noise in the power spectra. Note that in these two power spectra, 
the influence of the n7i-artefacts are still present even in the error 
weighted standard fit RM map. 

The two remaining power spectra in Fig. fallow us to investi- 
gate the influence of the nn- artefacts. The dashed-dotted line rep- 
resents the power spectra as calculated from a Pacman RM map 
allowing cr™ ax = 20° performing an error weighted fit. Again, no 
window weighting was applied to this calculation. One can clearly 
see that there is one order of magnitude difference between the er- 
ror weighted standard fit power spectra and the Pacman one at large 
k- values. Since we also allow higher noise levels cr™ ax for the po- 
larisation angles, this drop can only be explained by the removing 
of the nn- artefacts in the Pacman map. That there is still a lot of 
noise in the map on small real space scales (large fc's), which gov- 
erns the power spectra on these scales, can be seen from the dotted 
power spectrum, which was determined as above while applying a 
simple window weighting in the calculation. There is an additional 
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Figure 3. Various power spectra for the cluster Abell 2255 determined from 
different RM map making scenarios are shown. The solid line represents a 
power spectrum calculated from a standard fit RM map with c™ ax = 20° 
where no error weighting was applied. The dashed power spectrum was ob- 
tained from a standard fit RM map with c™ ax = 10° where error weight- 
ing was applied to the standard fit. The dashed-dotted power spectrum was 
determined from a Pacman RM map with cr™ ax =20° and error weight- 
ing was applied to the RM fits. The dotted line represents a power spectrum 
calculated from a Pacman RM map as above, but additionally a simple win- 
dow weighting was applied for the determination of the power spectra. For 
comparison, the power spectra of pure white noise is plotted as a straight 
dashed line. 

difference of about one order of magnitude between the two power 
spectra at large k- scales determined for the two Pacman scenarios. 

For comparison the power spectrum of pure white noise is 
plotted as straight dashed line in Fig. [3] It can be shown analyt- 
ically, that white noise as observed through an arbitrary window 
results in a spectrum of SB(k) oc k 2 . 

As another independent statistical test, we applied the gradi- 
ent vector alignment statistic V to the RM and (p° maps calculated 
by the Pacman and the standard fit algorithm. The gradients of RM 
and c p° were calculated using the scheme as described in footnote 
5 m Enfil in etalJ (Eo03). For the two RM maps shown in Fig.[Q we 
find a ratio of Ktdfit/V^acman = 20, which indicates a significant 
improvement mostly resulting from removing the nn- artefacts by 
the Pacman algorithm. The calculation of the normalised quantity 
V yielded V^tdfit = —0.75 and V p acman — -0.87. Smoothing the 
Pacman RM map slightly leads to a drastic decrease of the quan- 
tities V and V. This indicates that the statistic is still governed by 
small scale noise. This can be understood by looking at the RM 
maps of Abell 2255 in Fig. [I] The extreme RM values are situated 
at the margin of the source which is, however, also the noisiest re- 
gion of the source. Calculating the normalised quantity for the high 
signal-to-noise region yields Vpacman = —0.56. For comparison, 
the normalised quantity for the standard fit of the same high signal- 
to-noise region is V^tdfit = —0.94. 

3.2 Hydra North 

The polarised radio source Hydra A, which is in the centre of the 
Abell cluster 780, als o known as the Hydra A clu ster, is located at 
a redshift of 0.0538 Jde Vaucouleurs et aljfl99lh . The source Hy- 
dra A shows an extended, two-sided radio lob e. Detailed X-ray 
studies have been performed on this cluster (e.g. llkebe et all 1997: 
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Figure 4. A pixel-by-pixel comparison for the Hydra north lobe of 
i^Mgtdfit values obtained using the standard fit plotted on the y-axis versus 
the RMp acman values calculated by employing our new algorithm Pac- 
man. The parallel scattered points at ± 3000-4000 rad m~ 2 and ± 10 000 
rad m -2 are a result of spurious solutions to the n7t-ambiguity obtained 
by the standard fit algorithm. For the calculation of the RM maps, all five 
frequencies were used. 



Pere s et ailll998l : lDavid et ail feoOP. They revealed a strong cool- 
ing core in the cluster cen tre. The Faraday rotatio n structure was 
observed and analysed bv iTavlor & Perlevl {1993). They reported 
RM values ranging between —1 000 rad m -2 and +3 300 rad m -2 
in the north lobe and values down to —12 000 rad m -2 in the south 
lobe. 

Polarisation angle maps and their error maps for frequencies 
at 7815, 8 051, 8 165, 8 885 and 14 915 MHz resulting from ob- 
servation with the Very Large Array were kindly provided to us by 
Greg Taylo r. A detailed descrip tion of the radio data reduct ion can 
be found in lTavlor et all J199CT) and lTavlor & Perlevl Jl993h . 

In this section, we concentrate on the north lobe of Hydra A 
and discuss the south lobe separately in Sect. 13.31 The south lobe 
is more depolarised, leading to a lower signal-to-noise than in the 
north lobe. The north lobe is a good example of how to treat noise 
in the data, while the south lobe gives a good opportunity to discuss 
the limitations and the strength of our algorithm Pacman. 

Using the polarisation data for the four frequencies at around 
8 GHz, we calculated the standard fit RM map which is shown 
in the upper middle panel of Fig. [l] The maximal allowed er- 
ror in polarisation angle was chosen to be cr™ ax = 25° and 
RM max — 15 000 rad m~ 2 . A Pacman RM map is shown in the 
middle of the middle panel in Fig. [l] below the standard fit map. 
This map was calculated using the four frequencies at around 8 
GHz with a (Jk^gGHz — 30° and additionally where possible the 
fifth frequency at 15 GHz was used with a cr™=i 5GH z = 35°. 
RM max was chosen to be 15 000 rad m -2 . From a visual compar- 
ison of the two maps, one can conclude that the standard fit map 
looks more structured on smaller scales and less smooth than the 
Pacman map. This structure on small scales might be misinter- 
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Figure 5. Various power spectra calculated for the RM maps for the north 
lobe of Hydra A are shown. The dashed dotted line and the dotted line rep- 
resent the power spectra calculated from the standard fit RM map and the 
Pacman RM map, respectively, while no window weighting was applied. 
A threshold window weighting assuming a cr^ M = 75 rad m -2 yield the 
dashed line power spectrum for the standard fit map and the solid line for the 
Pacman map. The standard fit maps were calculated from the four frequen- 
cies at 8 GHz whereas the Pacman maps were determined using additionally 
the fifth frequency when possible. For both algorithms, cr™ ax = 35° was 
used. The influence of the small-scale pixel noise on the power spectra can 
clearly be seen in this figure at large k -values. 



preted as small scale structure of the RM producing magnetic field. 
Note that the difference in the RM maps is mainly due to the us- 
age of the fifth frequency for the Pacman map demonstrating that 
all information available should be used for the calculation of RM 
maps in order to avoid misinterpretation of the data. 

Since the four frequencies around 8 GHz are close together, a 
wrong solution of the rm-ambiguity would manifest itself by dif- 
ferences of about ARM = 10 000 rad m ~ 2 . Such jumps are not 
observed in the main patches of the standard fit RM map shown 
in Fig.[Q However, including the information contained in the po- 
larisation angle map of the 5th frequency at 15 GHz, which is de- 
sirable as explained above, one introduces the possibility of nn- 
ambiguities resulting in ARM = 3 000. ..4 000 rad m -2 . There- 
fore, a scatter plot of a pixel-by-pixel comparison between a stan- 
dard fit and a Pacman map, calculated both using the additional 
available information on the fifth frequency, is shown in Fig. [4] 
The parallel scattered points are spurious solutions found by the 
standard fit algorithm. One can clearly see that they develop at ± 
3000-4000 rad m -2 and less pronounced at ± 10 000 rad m -2 . 

As we will discuss below, the 15 GHz data set may have a 
lower signal-to-noise ratio than the 8 GHz data and thus, a lower 
error threshold cr™ ax will eliminate most of the wrong fits. How- 
ever, it is a particular strength of Pacman that it yields still reliable 
results even if choosing larger error thresholds. 

In order to study the influence of the noise on small scales, 
we calculated the power spectra from RM maps obtained using dif- 
ferent parameter sets for the Pacman and the standard fit algorithm. 
As for Abell 2255E, all calculations of power spectra were done for 
RM areas which would be covered by the Pacman fit if the same 
parameters were used. Again this leads to exclusion of pixels in the 
standard fit RM map which are not associated with the source. 

There is a clear depolarisation asymmetry of the two lobes 
of Hydra A observed as described by the Laing-Garrington effect 
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(Taylor & Perley|l993|;|Lane et al|2004|)- Therefore, we assume for 
the calculation of any power spectra for the Hydra source, that the 
source plane is tilted by an inclination angle of 45° where the north 
lobe points towards the observer and the south lobe away from the 
observer. 

For a first comparison, we calculated the power spectra for 
two RM maps similar to the one shown in Fig. [l] They were ob- 
tained using the four frequencies at 8 GHz for the standard fit al- 
gorithm and using additionally the fifth frequency when possible 
for the Pacman algorithm. The other parameters were chosen to 
be cr^ ax = 35° and RM max = 15 000 rad m~ 2 . For these two 
RM maps, we determined power spectra applying firstly no win- 
dow weighting at all and secondly threshold window weighting as 
described in the end of Sec. 12.31 We choose the threshold to be 
crf M = 75 rad m - 2 which represents the high signal-to-noise re- 
gion. The respective spectra are exhibited in Fig. [5] 

For the power spectra, nn- artefacts should play only a mi- 
nor role in the calculation of the power spectra as explained above. 
Therefore, any differences arising in these spectra are caused by the 
varying treatment of noise in the map making process or in the anal- 
ysis. Comparing the power spectra without window weighting, the 
one obtained from the Pacman RM map lies below the one from the 
standard fit RM map. The difference is of the order of half a magni- 
tude for large k's- small real space scales. The standard fit power 
spectrum seems to increase with larger k but the Pacman one seems 
to decrease over a greater range in k. The reason for this difference 
is the same one which was responsible for the different smoothness 
in the two RM maps shown in Fig. [T] namely that only four fre- 
quencies were used for the standard fit RM map. Even introducing 
a window weighting scheme, which down-weights noisy region, 
cannot account entirely for that difference as can be seen from the 
comparison of the window weighted power spectra. This is espe- 
cially true because the small-scale spatial structures are found even 
in the high signal-to-noise regions when using only four frequen- 
cies for the RM map calculation. This demonstrates again how im- 
portant it is to include all available information in the map making 
process. 

In order to study the influence of the maximal allowed mea- 
surement error in the polarisation angle cr™ ax , we calculated the 
power spectra for a series of Pacman RM maps obtained for cr™ ax 
ranging from 5° to 35°. The RM maps were derived using all 
five frequencies k m in= 5, not allowing any points to be consid- 
ered for which only four frequencies fulfil . < cr™ ax . Again an 
i?M max = 15 000 rad m -2 was used. The resulting power spectra 
using a threshold window weig hting (erf M = 50 rad m~ 2 ) are ex- 
hibited in Fig.[6| One can see that they are stable and do not differ 
substantially from each other even though the noise level increases 
slightly when increasing the cr™ ax . This demonstrates the robust- 
ness of our Pacman algorithm. 

For comparison, we plotted two more power spectra in Fig.[6| 
For the solid red line spectrum, no window weighting was applied 
to an RM map which was calculated as above having cr™ ax = 35°. 
One can clearly see that the spectrum is governed by noise on large 
/c-scales - small real space scales. Thus, some form of window 
weighting in the calculation of the power spectra seems to be nec- 
essary in order to suppress a large amount of noise. 

Another aspect arises in the noise treatment if we consider 
the power spectra represented by the blue stars which was calcu- 
lated using a threshold window weighting (<ro = 50 rad m -2 ). This 
power spectrum represents the analysis for an RM map obtained 
using only the four frequencies at 8 GHz while cr™ ax = 35°. It is 
striking that the noise level on large k- scales (on small real space 
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Figure 6. Power spectra for the north lobe of Hydra A calculated from a se- 
ries of Pacman RM maps are shown. The various coloured dashed and dot- 
ted lines show power spectra from RM maps calculated for cr™ ax ranging 
from 5° to 35° using all five frequencies fc m i n = 5. The spectra were de- 
termined applying threshold window weighting (erf M = 50 rad m -2 ). For 
comparison the solid red line power spectrum, which is the largest spectrum 
at k > 10 kpc -1 , is not window weighted for <j™ ax = 35°. The power 
spectrum represented by the blue stars was derived from an Pacman RM 
map obtained using the four frequencies at 8 GHz and a cr™ ax = 35°. 

scales) is lower than for the RM maps obtained for the five fre- 
quencies, one would expect it to be the other way round. Appar- 
ently the fifth frequency has a lot of weight in the determination of 
cr^ M (see Eq. (3) in paper I) since this frequency is almost twice 
as large as the other four. If the measurement errors of the polarisa- 
tion angles are underestimated this can lead to an underestimation 
of the uncertainty cr^ M in the final RM value due to Gaussian er- 
ror propagation. A more accurate error estimate cou ld be achieved 
by con sidering the error treatment as described by I Johnson et alJ 
( 1995). They state in their Appendix that errors are not uniform 
across synthesis images and propose a different error algorithm. 

One can test the hypothesis of underestimation of errors by 
performing a xl test as described in Sect. 12.21 The % 2 . map for 
the Pacman RM map of Hydra North derived by using a multi- 
frequency fit (see middle upper panel of Fig. [I) is exhibited in the 
left panel of Fig. [7] The red values are close to unity. However the 
yellow and blue regions in this map are values larger than unity in- 
dicating that the errors are underestimated in these regions. Note 
that in these regions there is also the fifth frequency used for the 
RM fit which might indicate that the errors for this frequency are 
underestimated by a larger factor than for the other four frequen- 
cies. 

We also applied the gradient vector product statistic V to the 
RM and (p° maps obtained by using the RM maps as shown in the 
middle panel of Fig.[I]and the respective ip° maps. The calculation 
yielded for the ratio V^tdfit /^pacman = 1.2. This corresponds to a 
slight decrease of correlated noise in the Pacman maps. The calcu- 
lation of the normalised quantity resulted in V^tdfit = —0.85 for 
the standard fit and Vpacman = —0.83 for the Pacman maps. This 
is expected since the low signal-to-noise areas of the north lobe of 
Hydra A will govern this gradient statistic and for both RM maps 
these areas look very similar. In order to visualise this, we calcu- 
lated a V map which is exhibited in Fig. The extreme negative 
values observed for the low signal-to-noise regions (red regions) in- 
dicate strong anti-correlation and thus, anti-correlated fluctuations 
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Figure 7. On the left hand side a xtij * s snown - The yellow and green regions have values larger than unity indicating underestimation of errors. Note that 
for these regions the RM fit used five frequencies whereas in the outer parts of the source only four frequencies were used. A V map for the north lobe of 
Hydra A is exhibited on the right panel. The extreme negative values (red and black) correspond to strong anti-correlation of RM and tp° fluctuations, which 
are produced by noise in the observed maps. However, the green and blue coloured regions have values around zero and are therefore high signal-to-noise 
regions (compare with the <r RM map in the middle panel of Fig.[TJ. 



of RM and ip° in these regions. For the high signal-to-noise regions 
(blue and green) moderate values varying around zero are observed. 
Note the striking morphological similarity of the RM error («7 RM ) 
map in the lower middle panel of Fig.[l]and of the V map as shown 
in Fig. These two independent approaches to measure the RM 
maps accuracy give a basically identical picture. These approaches 
are complementary since they use independent information. The 
RM error map is calculated solely from <jk maps, which are based 
on the absolute polarisation errors, whereas the V map is based 
solely on the RM derived from the polarisation angle ip(k) maps. 



However, we note that the resul ts for magnetic field s trengths 
and correlation lengths presented bv lVost & E nfilin (2003) ar e not 
changed substantially. This is owing to the fact that Vogt & Enfilin 
(2003) considered for the calculation of the magnetic field proper- 
ties only the power on scales which were larger than the resolution 
element (beam). Therefore the small scale noise was not consid- 
ered for this calculation and thus, no change in the results given 
is expected. The central magnetic field strength for this cluster as 
derived during the course of this work is about 10 \iG and the field 
correlatio n length is about 1 kp c being consistent with results pre- 
sented in lVogt & Enfiml fcOOfl . 



3.3 The Quest for Hydra South 

As already mentioned, the southern part of the Hydra source has to 
be addressed separately, as it differs from the northern part in many 
respects. One of them is that the depolarisation is higher in the 
south lobe which especially complicates the analysis of this part. 

An RM map obtained employing the standard fit algorithm us- 
ing the four frequencies at 8 GHz is shown in the upper right panel 
ofFig.[T]This RM map exhibits many jumps in the RM distribution 
which seems to be split in lots of small patches having similar RMs. 
Furthermore, RMs of —12 000 rad m -2 are detected (indicated as 
red regions in the map). If one includes the fifth frequency at 15 
GHz in the standard fit algorithm, the appearance of the RM map 
does not change significantly, and most importantly, the extreme 
values of about —12 000 rad m -2 do not vanish. 

The application of the Pacman algorithm with conservative 
settings for the parameters cr™ ax , cr^ ax and g, leads to a split- 
ting of the RM distribution into many small, spatially disconnected 
patches. Such a map does look like a standard fit RM map and there 
are still the RM jumps and the extreme RM values present. How- 
ever, if one lowers the restrictions for the construction of patches, 
the Pacman algorithm starts to connect patches to the patch with the 
best quality data available in the south lobe using its information 
on the global solution of the nn- ambiguity. An RM map obtained 
pushing Pacman to such limits is exhibited in the middle right panel 
of Fig. [I] Note that the best quality area is covered by the bright 
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Figure 9. Various power spectra for the RM maps of the south lobe of the 
Hydra source are compared with the power spectrum of the north lobe cal- 
culated from the Pacman RM map (fc m i n = 5, cr™ ax = 35°) employing a 
threshold weighting with cr^ M = 50 rad m -2 represented as filled circles. 



by comparing the threshold window weighted power spectra (dot- 
ted line for the power spectra of the Pacman RM map and dashed 
dotted line for the power spectra from the standard fit RM map) in 
Fig- El Since the choice of the window weighting scheme seems to 
have also an influence on the result the situation is still inconclu- 
sive. 

An indication for the right solution of the n7r-ambiguity prob- 
lem may be found in the application of our gradient vector product 
statistic V to the Pacman maps and the standard fit maps of the 
south lobe (as shown in Fig. [TJ. The calculations yielded a ratio 
V^tdfit/Vpacman = 12. This is a difference of one order of magni- 
tude and represents a substantial decrease in correlated noise in the 
Pacman maps. The calculation of the normalised quantity yielded 
V'stdfit = —0.69 for the standard fit maps and Vpacman = —0.47 
for the Pacman maps. This result strongly indicates that the Pacman 
map should be the preferred one. 

The final answer to the question about the right solution of 
the RM distribution problem for the south lobe of Hydra A has to 
be postponed until observations of even higher sensitivity, higher 
spatial resolutions and preferably also at different frequencies are 
available. 



blue regions and a zoom-in for this region is also exhibited in the 
small box to the lower left of the source. In the lowest right panel 
of Fig.[l] a a RM -map is shown indicating the high quality regions 
by the red colour. 

The RM distribution in this Pacman map exhibits clearly 
fewer jumps than the standard fit map although the jumps do not 
vanish entirely. Another feature that almost vanishes in the Pacman 
map are the extreme RMs of about — 12 000 rad m -2 . The RM dis- 
tribution of the Pacman map seems to be smoother than the one of 
the standard fit map. However, if the individual RM fits for points, 
which deviate in their RM values depending on the algorithm used, 
are compared in a ip(k)- A 2 -diagram, the standard fit seems to be 
the one which would have to be preferred since it does fit better to 
the data at hand. Pacman has some resistance to pick this smallest 
X 2 solution if it does not make sense in the context of neighbour- 
ing pixel information. As an example, we plotted individual RM 
fits which are observed across the RM jumps in Fig. [8| However 
as shown above, one has to consider that the measurement errors 
might be underestimated which gives much tighter constraints on 
the fit than it would otherwise be. 

In order to investigate these maps further, we also calculated 
the threshold weig hted (cr? M = 50 rad m~ 2 ) and simple weighted 
power spectra from the Pacman and the standard fit RM maps of 
the south lobe and compared them to the one from the north lobe. 
The RM maps used for the comparison were all calculated employ- 
ing cr™ ax = 35° and k m in = / = 5. The various power spectra are 
shown in Fig.[5| One can clearly see that the power spectra calcu- 
lated for the south lobe lie well above the power spectra from the 
north lobe, which is represented as filled circles. 

Concentrating on the comparison of the simple window 
weighted power spectra of the south lobe in Fig.[9|(solid line from 
Pacman RM map and dashed line from standard fit RM map) to the 
one of the north lobe, one finds that the power spectrum calculated 
from the Pacman RM map is closer to the one from the north lobe 
than the power spectrum derived from the standard fit RM map. 
This indicates that Pacman RM map might be the right solution to 
the RM determination problem for this part of the source. However, 
this difference vanishes if a different window weighting scheme is 
applied to the calculation of the power spectrum. This can be seen 



4 CONCLUSIONS & LESSONS 

We demonstrated the robustness of our Pacman algorithm which 
is especially useful for the calculation of RM and cp° maps of ex- 
tended radio sources. It dramatically reduces nTi-artefacts in noisy 
regions and makes the unambiguous determination of RMs in these 
regions possible. Any statistical analysis of the RM maps will profit 
from this improvement. 

In the course of this work the observational data taught us the 
following lessons: 

(i) It is important to use all available information obtained by the 
observation. It seems especially desirable to have a good frequency 
coverage. The result is a smoother, less noisy map. 

(ii) For the individual least- squares fits of the RM, it is advisable 
to use an error weighting scheme in order to reduce the noise. 

(iii) Global algorithms are preferable, if reliable RM values are 
needed from low signal-to-noise regions at the edge of the source. 

(iv) Sometimes it is necessary to test many values of the parame- 
ters RM max , cr™ ax and <7m ax , which govern the Pacman algorithm, 
in order to investigate the influence on the resulting RM maps. 

(v) One has to keep in mind that whichever RM map seems to 
be most believable, it might still contain artefacts. One should do 
a careful analysis by looking on individual RM pixel fits but one 
also should consider the global RM distribution. We presented the 
gradient vector product statistic V as a useful tool to estimate the 
level of reduction of cross-correlated noise in the calculated RM 
and ip° maps. 

(vi) For the calculation of power spectra it is always preferable 
to use a window weighting scheme which has to be carefully se- 
lected. 

(vii) Finally, under or overestimation of the measurement errors 
of polarisation angles (which will propagate through all calcula- 
tions due to error propagation) could influence the final results. We 
performed a reduced xl test m order to test if the errors are over or 
under estimated. 

We conclude that the calculation of RM maps is a very dif- 
ficult task requiring a critical view to the data and a careful noise 
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Figure 8. An example for individual RM fits across an RM jump observed in the RM maps of the south lobe of Hydra A. The solid line indicates the standard 
fit solution and the dashed one represents the Pacman solution. The error bars in these plots indicate the polarisation errors which were multiplied by 3. Note 
that this is an extreme example, where Pacman chooses a solution to the individual fits in noisy regions which has not necessarily the minimal x 2 value by 
construction. However, the global statistical tests applied to the whole RM map of Hydra South indicate that the Pacman map has less artefacts in comparison 
to the standard fit RM map. We note, that for this particular local example it is not trivial to decide which algorithm is giving the right solution. 



analysis. However, we are confident that Pacman offers a good op- 
portunity to calculate reliable RM maps and to understand the data 
and its limitations better. Furthermore, such improved RM maps 
and the according error analysis for Hydra North and Abell 2255 
were presented in this paper. 
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